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Abstract. Boundary effects in quantum mechanics are examined by considering a 
partition wall inserted at the centre of a harmonic oscillator system. We put an equal 
number of particles on both sides of the impenetrable wall keeping the system under 
finite temperatures. When the wall admits distinct boundary conditions on the two 
sides, then a net force is induced on the wall. We study the temperature behaviour 
of the induced force both analytically and numerically under the combination of the 
Dirichlet and the Neumann conditions, and determine its scaling property for two 
statistical cases of the particles: fermions and bosons. We find that the force has 
a nonvanishing limit at zero temperature T — and exhibits scalings characteristic 
to the statistics of the particles. We also see that for higher temperatures the force 
decreases according to 1/y/T, in sharp contrast to the case of the infinite potential 
well where it diverges according to VT. The results suggest that, if such a nontrivial 
partition wall can be realized, it may be used as a probe to examine the profile of the 
potentials and the statistics of the particles involved. 
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1. Introduction 

Quantum systems are often delineated by modeling their classical counterparts — in 
fact, it is a standard practice that we define a system in quantum mechanics through the 
procedure called 'quantization', which amounts to replacing functions of phase space in 
classical mechanics by appropriate operators based on commutation relations. However, 
this quantization procedure does not necessarily provide a unique quantum system 
to a given classical system, with the familiar example being the ordering ambiguity 
of operators. Nontrivial topology of the classical configuration space furnishes an 
additional ambiguity in the quantum system, which is exemplified by a particle moving 
on a circle where a multiple of windings are allowed for transitions. In particle physics, 
the same topological effect is known to be responsible for the infinite vacua structure 
which causes the strong CP violation (see, e.g., pQ). 

Another source of ambiguity, which is less recognized and yet physically more 
tangible than the aforementioned ones, lies in the choice of boundary conditions imposed 
on quantum states [2j [3l 0]. A prime example of this may be found in a particle system 
with an impenetrable wall, where one has a variety of boundary conditions to choose, 
say, from Dirichlet to Neumann or anything in between. These boundary conditions 
represent different physical properties of the wall which are missing in the classical 
description, and can lead to novel effects in the quantum system j5] . A further example is 
provided by a singular point on a line, which serves as a source for a number of interesting 
properties, including supersymmetry and Berry phase [Tl, [8] . The importance of 
boundary conditions in quantum mechanics becomes evident if we recall the rudimentary 
fact that a different choice of boundary conditions yields a different spectrum. The fact 
that the physical properties of quantum dots depend heavily on the boundary conditions 
will also be worth mentioning. 

In order to provide a simple setup where the boundary effects can be observed 
directly, in our previous work j9] we presented a case study of the induced pressure, or 
statistical quantum force, which emerges on an impenetrable partition wall inserted at 
the centre of an infinite potential well, when the wall is assumed to realize a nontrivial 
set of boundary conditions: the Dirichlet condition on one side and the Neumann on the 
other. We will not delve into here how such a partition can actually be manufactured as 
a device, but only mention that walls admitting generic boundary conditions including 
Dirichlet and Neumann can be realized by a combination of square well potentials in 
the vanishing limit of their widths [5]. When two such walls, one with Dirichlet and 
the other with Neumann, are 'glued together' within a narrow distance, then it will 
serve as our partition effectively under a scale significantly larger than the distance. In 
more formal terms, our partition is a special example of the general {U(2) family of) 
point singularities allowed quantum mechanically on a line, whose realizations by scaled 
families of regular potentials has been studied extensively [I]. 

Once the partition wall is realized and placed in the well separating the same 
number N of particles on its two sides, under the finite temperature T we expect that 
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Figure 1. The eigenfunctions and eigenvalues under the harmonic potential V(x) — 
^-x 2 with a partition at the centre. If the partition admits the Dirichlet (ip = 0) and 
the Neumann (?/>' = 0) boundary conditions on the left and on the right, respectively, 
the eigenstates (the lowest two are shown in both half lines) possess different energy 
levels. When the same number N of particles are introduced in each of the half lines, 
these level differences give rise to a net force on the partition. 



the wall is pushed from the two sides by the particles in motion. Now, the point is that 
the different boundary conditions imposed at the partition cause different energy levels 
and, accordingly, different statistical distributions of particles in the levels between the 
two half wells, yielding a net force AF on the wall. In [9J we investigated the behaviour 
of the force AF as a function of T, and examined how it scales with particle number N 
for the two kinds of particle statistics, the Bose-Einstein statistics and the Fermi-Dirac 
statistics. There we have found that the force AF has a finite limit at zero temperature 
T = which scales as N for bosons and as iV 2 for fermions, and that it has a minimum 
before it diverges as VT for the high temperature limit T — > oo. 

This raised a natural question if these results are specific to the potential used in the 
analysis, and if so how. In the present paper, we attempt to answer this by considering 
the same partition wall placed in a different potential, namely, the harmonic oscillator 
potential (see Figured]). One of the reasons for the harmonic potential is that, unlike 
the infinite well potential, it stretches infinitely for higher energies and shares a feature 
with potentials which are often used to describe actual physical systems. Another reason 
is that, on account of its technical simplicity which we also exploit here, the harmonic 
potential itself is widely used in various physical contexts including confinement of 
particles in a narrow region. Again, as illustrated in Figure [2j we consider both bosonic 
and fermionic cases for particle statistics, and the number N of particles is regarded as 
large but not macroscopically large (to ensure its possible relevance to nano devises). 
We derive analytic approximate formulae for the force in the low, medium and high- 
temperature regimes separately, which are found to be reasonably good to reproduce 
the numerical results obtained for N = 100. Note that, in actual realizations, our 
one dimensional system can be regarded as a model of an axis in three dimensions 
perpendicular to the surface of the partition which attracts the particles by the harmonic 
potential. 

Our analysis then shows that, in contrast to the infinite well case, the force AF 
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Figure 2. Illustration of particle distributions over the levels in the two half lines 
at a low temperature. The force related to each level differs on the two sides of the 
partition and induces a non-vanishing net force on the partition, which is dependent 
on temperature as well as on the particle statistics. 



on the wall decreases as 1/VT as the temperature increases in the high temperature 
regime, and eventually vanishes in the limit T — > oo. This is the case for both bosons 
and fermions, and the force is of the order of N. For the low temperature regime, the 
zero temperature limit of the force is of the order of iV for bosons as in the potential 
well case, but for fermions it is of the order of y/~N in contrast to iV 2 obtained in the 
potential well. We also find that, unlike the infinite well case, no minimal point of the 
net force is found in the medium-temperature regime irrespective of the statistics of 
the particles. Implications of these results, combined with the previous ones, will be 
discussed in the text. 

This paper is organized as follows. In section [2] we define the model and 
provide our scheme of analytical approximation for the induced force at arbitrary finite 
temperatures. In section [31 we present our analysis in detail for the high temperature 
regime as well as the numerical results obtained. Section 0] is devoted to the analysis of 
the low temperature regime, where we employ independent approaches for the fermionic 
and bosonic cases. Section [5] discusses the medium-temperature regime where we also 
seek to interpolate the other two temperature regimes. Finally, we present our conclusion 
and discussions in section [6j 

2. Definition and basic properties of the system 

In this section, we first define the system mentioned earlier and introduce some 
notations convenient for the description of the system. With these, we present our 
basis of evaluating the induced force on the wall, and provide our scheme of analytical 
approximation used later. 

2.1. The eigenstates and distributions at finite temperatures 

To begin with, we recall the basic result of a harmonic oscillator system. For a particle 
moving on a line under the harmonic potential V(x) = t^rf-x 2 , the normalized energy 
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eigenf unctions 4>i( x ) an d the corresponding energy eigenvalues £ obeying the equations, 

. . , . h 2 d 2 muj 2 9 , . 

**<«)=**<«), H=--- i+ — X \ (1) 

are given by 

^H^i^^Qe"^, £ l= i (/ + £), / = 0,1,2,.... (2) 



Here, Hj(^) are the Hermite functions [TO] and we have introduced the length scale 
a = \ — and the energy scale £ = 2huo = (the unusual factor 2 is for our 

later convenience). The Hermite functions are even functions for even I = 2k and odd 
functions for odd I — 2k + 1, and fulfill 

H 2fc (0) = (-1)*^, H' 2fc+1 (0) = 2(2fc + l)H 2fe (0). (3) 

Now let us divide the line into two by inserting an infinitely thin wall at x = and 
consider its consequence in one of the half lines, say, the positive one. If the boundary 
condition on the wave functions at the wall is the Dirichlet condition, then for the 
particle confined in the positive half line, we have only the odd eigenfunctions, that is, 

<4 Dir) = V2&*+i, 4 Dir) = W, A; = 0,1,..., (4) 

where the factor v2 is required for the normalization in the half line. On the other hand, 
if it is the Neumann boundary condition, then we have only the even eigenfunctions, 

4 Neu) = V20 2 „ Et u) =£2 k , fc = 0,l,... (5) 

One may wonder if there exist eigenfunctions other than these fll]) and (jSJ) in our system, 
since after all our partition wall breaks the parity symmetry under x —x forcing us to 
work in two half lines separately, and there seems no reason to consider only the standard 
eigenfunctions of the harmonic oscillator defined on the whole line with definite parities. 
However, this cannot be the case, because if ip+(x) is any eigenfunction defined on the 
positive half line obeying, say, the Dirichlet condition at x — 0, then one can extend it to 
the whole line by setting ip{x) = ip+(x) for x > and ip(x) = —ip + (—x) for x < 0. The 
resultant function ip(x) is regular at x = and satisfies the equation ([1]) on the whole 
line and hence provides an eigenfunction with odd parity of the harmonic oscillator, 
implying that it should be one of the standard ones <p2k+i- The case of the solutions 
with the Neumann condition can be argued similarly. 

To proceed, consider a statistical system consisting of iV mutually noninteracting 
identical particles in the harmonic potential in the half line under finite temperature T. 
For convenience, we introduce the dimensionless quantities, 

t = k B T/£, b = 1/t, (6) 
with k-Q being the Boltzmann constant, and write the energy eigenvalues as 

E k = £e k , with e k = (k + a), k = 0, 1, . . . , (7) 
where the constant a is given, for the Dirichlet case and the Neumann case, by 

a (Dir) = 3 ^(Neu) = \ (g) 

4 4 
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respectively. The statistical distributions of the particles then read 

1 1 f 1 bosons /nS 

Nk = = n — < (9) 

e a+be k _ Tj e ( a +ba)+bk _ [ -1 fermions 

At a given temperature, the chemical constant a is determined uniquely from 
the constraint N = Y^k=o^k- This implies that the combination a + ha is uniquely 
determined by the total number constraint, from which we learn that the distributions 
are actually the same for both of the Dirichlet and the Neumann boundary conditions, 
that is, N^ 1 ^ = Nj^ cn ^ for all k. This allows us to introduce 



(Dir) + ^(Dir) = ^(Neu) + ^(Neu^ ^ . = ^(Dir) = ^(Ncu) ^ 



a := a 

to obtain the simpler expression, 

N ^ = ^—^ fc = o,i,..., (11) 

commonly used for the two boundary conditions. We also find from l/iV fc = e a+bk — n 
and b > that 1/JV < l/Nt < 1/N 2 < ■ ■ ., or JV > N x > N 2 > . . .. Note that for 
bosons 77 = 1, the positivity of distributions iVi, > for all k implies a > at any 
temperature (where the case k = gives the strongest condition). For fermions 77 = — 1, 
no such restriction emerges and a can take any value in (—00, 00). 

At this point, we mention that the sum over the levels k admits an exact 
resummation valid for a > 0, 



00 1 00 (a+6fc) 00 a bk 

N = V - = V — = 7iY — 

^-^ oa+fefc r, ^— ' 1 tip - (a+bk) 1 1 ^.p—a^—bk 

k=0 C fc=0 L fc=0 1 e 



OO OO ,O0 ,O0 OO ( flQ &\ 

= V E E (^e- 5 e- 6fe ) = , E fa"*) E ^ = V E ^ ' (12) 

fc=0 2=1 2=1 fc=0 1=1 e 

where we have used n^ 1 = i]. This resummation formula will be useful later. 



2.2. The force difference 

Let us now suppose that the wall inserted at x = in the harmonic potential imposes 
the Dirichlet boundary condition on the left (negative) side and the Neumann boundary 
condition on the right (positive) side. Due to the difference in the energy levels developed 
in the two sides of the wall, one expects that a net force, or statistical pressure, will 
emerge on the wall as a purely quantum effect deriving from the boundary conditions 
(see Figure [2]). Our aim is to evaluate this induced net force as a function of (rescaled) 
temperature t. 

Before we proceed, we recall the fact that for the case of the infinite potential well 
[9], the force acting on the wall from each side of the half lines proves to be essentially 
the same as the one giving the average energy, E = Y^k=o -^kEk- From this, the net 
force is obtained by the difference AE = £'( Dir ) — E^ Nen \ For the harmonic oscillator, 
the difference in the average energy reads 

oo c oo Arc 

AE = Y,N k (Ef " 4 NCU) ) = f E *k K ir) - - (Ncu) ) = ^, (13) 

fc=0 1 k=0 1 
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which is temperature- independent. 

However, for the harmonic oscillator the average energy is no longer the same as 
the force. To see this, let us consider the contribution F for the total force coming from 
one single level E. Under a shift 5x of the wall from the origin, the energy level will 
also change by 5E, and from this the force is found by 



which is the difference of the forces between the two sides of the wall arising from the 
two corresponding energy levels specified by the same integer k. Unfortunately, unlike 
the infinite well case, we do not have analytical solutions for the harmonic oscillator 
when the partition is displaced from the centre, and we are compelled to resort to some 
approximation scheme to evaluate the force (I14p . 

At this point we recall that, mathematically speaking, the Hamiltonian operator H 
of our system has the infinity x = ±00 as a limit-point singularity whereas the position 
of the wall is a regular endpoint, meaning that for any real eigenvalue E, there exists 
only one normalizable eigenfunction up to a phase factor. In more concrete terms, given 
an arbitrary E we have two independent solutions for the differential equation (pQ) but 
requirement of normalizability determines a particular linear combination of the two 
as a possible candidate for an eigenfunction. It qualifies as a true eigenfunction when 
the boundary condition at the partition is further met, which is attained by tuning E 
to be one of the particular set of real numbers which form the energy spectrum of the 
system. This heuristic picture of approaching eigenfunctions and eigenvalues suggests 
that, if the shift 5x of the wall is sufficiently small, for a fixed (Dirichlet or Neumann) 
boundary condition the difference in the eigenfunctions should be small in the L? sense 
that their scalar product tends to 1 as 5x — > 0, with the perturbed eigenvalue E also 
being close to the unperturbed one E. Similarly, we expect that for a small variation 
Sx the difference in the boundary values of the two wave functions or their derivatives 
- the former is nonvanishing for the Neumann case while the latter is nonvanishing for 
the Dirichlet case - remains small and in the same order Q(8x) at most. 

Now, we consider the identity valid for any two real and normalized eigenfunctions, 



where ( • , • ) + denotes the scalar product on the positive half line, and the prime indicates 
the derivative with respect to x, e.g., ip' = -^ip. Specifically, for the Dirichlet case, we 
choose in (115j) an unperturbed Dirichlet eigenvalue E^ ir ^ and its eigenfunction f^ 1 ^ for 
E and (p E , and the perturbed eigenvalue and eigenfunctions caused by the shift in the 
wall for E and (fE, respectively. Based on our observations on the perturbed quantities, 
we find that the formula ( |T5|) in the leading order of 5x or of 5E = E — E^ 1 ^ yields 




(E-E) ((p E ,(p E ) + = ((p E ,H(p E ) + - (H(p M ,(p E ) 



h 2 r 

= — ^(0)^(0) -^(0)^(0) 
2m 1 



(15) 



6E = E - 4 Dir) 



2m 



^r j/ (o)^(o), 



(16) 
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In parallel, our assumption ensures that 

^(0)^f r) '(0), 
and that the Dirichlet condition is satisfied at the shifted wall x = 5x, 

= <p E (5x) « <^(0) + (p E '(0) ■ 5x, 
from which we have 



'5E\ iBiT) 



2m 



^ Dir)/ (o) 



Analogously, for the Neumann case, we have 

h 2 



or 



5E = E-E 



(Neu) 



^r(0)9.'(0), 



SE 



Since (p E is an eigenfunction, we find 

2m£ 



2m£ 



(Neu) 



M'(o) = — ^ — ^- 

The Neumann condition is satisfied at x = <5x if 
= <p E '(6x) « ^(0) + <^/'(0) • fe, 
from which we obtain 

fdE^ ■* ( Ncu ) r (Ncu)/ nX l 2 

[■n) t xE[ l« (0) 1 • 

Combining (120]) and (1251) together with 



(Neu) / n N 
¥4 (0)- 



and 



(17) 

(18) 
(19) 

(20) 

(21) 
(22) 

(23) 
(24) 

(25) 

one can evaluate the 



contribution for the force difference coming from the fc-th level as 

-2 



AF fc = 

2m 

h 2 



^f r)/ (o) 

2 



E 



(Neu) 



fk (o) 

H 2fe+1 (0) 



1 



2ma? y/^2 2k+1 (2k+ 1)! 
£ (2k)\ 



£ k 



4/ v^Fa 2 2fc (2£;)! 



[H 2fc (0)f 



" 2v^Fa2 2fe (fc!) 2 ' 

The final expression shows that AFj. has a nontrivial /c-dependence in contrast to the 
infinite well case jH] where it is simply proportional to k. 

In what follows, for brevity we shall use the dimensionless force difference, 

AF fc (2A;)! 



A/ fc := 2^. 



ixa- 



2 2k (k\) 



(26) 
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which has the first few values, 

A/ = l, A/i = A/ 2 = | A/ 3 = ^- (27) 

Note that, in general, Af k > Af k+1 (since Af k+1 /Af k = (2k + 1) /(2k + 2) < 1) and, 
consequently, Af k < 1. To see how it behaves for large k, we may use the Stirling 
formula \U\ 



— ) \f2nri 



to obtain 



(2*Q! 
2 2k (k\) 



'nk 



1 



12n 
1 



n- 



8fc 



+ 



1 

A-2 



(28) 



(29) 



which provides a sufficiently good approximation already at k — 1, and improves quickly 
for larger One can also observe that the 'rearranged' approximation 

1 



A/ fc « 

\A (* + 1/4) 

is similarly good as ( 1291 . and can be used for = as well (see Table 1). 



(30) 





A/ fc 


i 

V 7rfe 


fl - 


i 

0r(fc+l/4) 





1 






1.12838 


1 


0.5 


0.56419 


0.49367 


0.50463 


2 


0.375 


0.39894 


0.37401 


0.37613 


3 


0.3125 


0.32574 


0.31216 


0.31296 



Table 1: Exact and approximate values for Af k , for k = 0, 1, 2, 3. 



The total force difference (the net force) is then given by 
00 00 (2kV 1 

y =S w =S?(^ (31) 

Now, using the resummation formula analogous to ( 1T2l together with the fact that a 
sum of the form 

is actually the Taylor expansion of 1/^/1 — q for |g| < 1, we finally arrive at the 
convenient analytical expression of the net force: 



(33) 
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We mention that the inequality Af k < 1 noted earlier implies 

oo oo 

Af = J2 N k Af k <Y,N k = N. (34) 

k=0 k=0 

Clearly, this inequality fl34|) is expected to be close to the equality for low temperatures 
but it will become loose as temperature increases where the energy levels of higher k are 
excited with increasing probability. This way we can expect that the net force decreases 
as temperature increases. One can also obtain an improved inequality utilizing the closer 
details A/ = 1 and Af k < 1/2 for k > 0, that is, 

oo oo oo 

Af = N k Af k = N - N k (l - Af k ) <N-Y,N k -l/2 = (N + N )/2. (35) 

k=0 k=l k=l 

We have furnished analytical approximations for the net force which form our basis 
for studying the temperature dependence of the force on the partition wall. We stress 
that, unlike the case of the potential well, we further need to improve the approximations 
to obtain results comparable to numerical computations. This is required by the 
nontrivial level dependence fl3T|) of the force in the harmonic case, and below we shall 
establish independent approaches to deal with infinite sums which are appropriate for 
three different temperature (high, low and medium) regimes. 



3. The high-temperature regime 

We now analyze the behavior of the net force when the temperature is sufficiently high 
t 3> 1 or b <C 1 (recall that t = 1/6 is the rescaled dimensionless temperature in (jSJ)), 
where we are allowed to take a > 0. One can then approximate the formula (fT2l as 

OO ( /Ig - "^ OO ( f)Q~ 

^ = ^E}^-^E i V ± - (36) 

1=1 1 e 1=1 Ul 

This approximation is certainly good for the terms I t, but if we assume the condition 
e~ a 1 (which is stronger than a > 0) for which higher-/ terms in the sum ( l36l) are 
suppressed, then we can rewrite ( )36|) and perform the summation in closed form as 

t l 1 — ne a 



i=i 



which implies 



a , N 7] /N\ 2 1 /NY ( (N s 



We can see from this outcome that for sufficiently high temperatures t ^> N our 
approximation is indeed consistent with e~ a <ti 1, suggesting that a valid approximation 
for a(t) can be obtained at least for t with t 3> N. Based on this observation, we shall 
define the high-temperature regime of the system by the condition t ^> N. Notably, 
here a depends on t and iV only through the combination t/N. 
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Figure 3. High temperature approximations for Af, for fermions (left) and bosons 
(right), N = 100. Thick solid line: the numerical result. Thin solid line: the expansion 
([40]) truncated to one term. Dashed line: truncation to two terms. Dash-dotted line: 
truncation to three terms. Dotted line: {43]). 



Similarly, we can also approximate the net force (1551) as 

00 (rie~ & ) 1 00 (rie~ & ) 1 _ / + \ 3 00 (rie~ 6! ) 1 

^"'gfebr-'gV-^^^) gV- (39) 



Inserting the expansion (1551) . we find 
1— fN 




(40) 

with the coefficients 

\/2-l 1-3^ + 2^ 11 + 7^2-12^3 
ci = ^— , c 2 = , c 3 = . (41) 

Their numerical values are estimated as 

ci « 0.207, c 2 w 0.0369, c 3 w 0.00479. (42) 

The analytical results, together with those improved below, are depicted in Figure [5] 
along with the numerical ones. There we find that, for both of the fermionic and 
the bosonic cases, the net force is a monotonically decreasing function in the high- 
temperature regime and vanishes in the limit t —* 00. In fact, the expression ( 1401) shows 
that it decreases according to 1/y/i asymptotically as t — > 00 with the common order 
N. This outcome is in sharp contrast to the infinite well case [9] where the force A/(t) 
diverges according to \/t for both fermions and bosons (see Figure H]) . 

It is worth noting that this high-temperature expansion can be made applicable even 
for medium temperatures if we modify it slightly by adopting the Pade approximant 
form [12] to regularize its diverging behavior for t — > 0. Explicitly, we may take 

/n\s (1)5 

Af w y/N > « r A f ~ ^ — yNJ r (43) 

x '/ 2 t n 'I 2 
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Figure 4. Numerical results of the net force Af in the infinite well case; for fermions 
(left) and bosons (right), at N — 100. The net force diverges as t — > oo according to 
y/i for both fermions and bosons [9]. 

both of which admit the expansion (BO]) but with slightly modified coefficients, 

ci « 0.207, c 2 « 0.0429, c 3 « 0.00888. (44) 

We note that the regularized formula reproduces the analytical approximation (1401) 
quite well up to the third term (the fourth also not being very different). However, 
for bosons the attempt to regularize at t — fails because (j4"5j) diverges at a certain 
positive medium temperature value, i.e., 4 = v ^~ 1 , and this also ruins improvement in 
precision. For fermions, in contrast, the approximation is indeed valid up to fairly low 
temperatures. It is also important to observe from the result (140"]) that, similarly to a, 
the net force Af/\/~N also depends on t and iV only through the combination t/N in 
the high-temperature regime. 

In passing, we provide a technical remark that, if one wishes to improve the 
approximations (136]) and (1391) to curb the error caused by the replacement 1 — e~ bl = 
bl + 0((bl) 2 ) « bl used above, one may instead use — e~% = bl + 0((bl) 3 ) ~ 6/ (which 
is more precise by one order) to obtain 

00 (rje~ & y 00 (rie~ &+ ^ 1 00 f^e"")' 00 (r/e^+ty 

S^^'S « ' Stt^'S vm ' (45) 

respectively. Although the acquired improvements in a and in Af are of the order of 
1/t and hence insignificant for iV — * oo, they may become significant for iV ~ 100. 

4. The low-temperature regime 

Next we turn our attention to the net force when the temperature is sufficiently low. 
In this regime we need to develop our approximation depending the statistics of the 
particles, and below we present our arguments for fermions and bosons, separately. 
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4-1. Fermions at low temperature 

Note first that, for the fermionic case 77 = —1, the particle distribution at t = becomes 
N — Nx — ■ ■ ■ — N N _, = 1, N N = N N+1 = ■ ■ ■ = 0, (46) 
for which a is, so to say, —00. The net force is then 



Af(0)- VAf V < 2fc > ! - 2N ( 2N V- 

n > h 22k ( kl ) 2 22Nm2 



1 "8iV +0 U 



(47) 



where the exact result is, at the last step, approximated for large N using the Stirling 
formula. Observe that the limiting value A/(0) is of the order \[N for large N. 

For a slightly higher temperature t > 0, we may assume that the distribution differs 
from (0SD only at N N ^ and N N . Then T,T= N k = N implies JVjv_i + A^ = 1, which 
can be exploited as 

1 1 e a+b(N-l) 

N N = l-N N ^ = l- ■ - = -g— ^ ■ - (48) 



^a+bN 

which is equivalent to 



e ^u„ + l "' JV x -'JV-i - e a+fe(V-l) + ]^ e a+6(JV-l) _|_ 1 



fi = -&(tf-i), (49) 

in this approximation. 

For the net force, we can use the same approximate distribution, but we may further 
suppose that (l4"9l) is a good approximation even when more than two levels — let their 
number be denoted by 2 J, still assuming 2 J <C 2N — get nontrivially occupied. Then 
the net force Af will differ from the zero temperature value A/(0) as 

00 JV-l N—l 00 

Af - A/(0) = ]T A/,iV fc - X) A /* = " E A ^ t 1 - ^) + E A /^ 

fc=0 fc=0 fc=0 fe=JV 

N-l N+J-l 

~- X A/jfc(l-iV fc )+ E A/ fc A^ fc 



k=N-J k=N 

1 



E a /at-j (1 - ^4 - ) + E A ^-l +J - 

i=i V e b U J J + 1/ i= i 



.KM) 



This can be further approximated by using ( 1301) as 



(50) 



J Af Af J -Z 3 (1 - 2j) 

A/ - A/(0) = E A/ ^,7 Af "~' « E , (51) 

j= i e V 3 2^ + 1 3=1 e y ] 2^ + 1 



from which we obtain 

Af « A/(0) 1 1 



6 1 36 

+ 1 e 2 + 1 



2^^ 

on account of j < J <C N. As an expansion in terms of e~^, (152]) turns into 



(52) 



i) 



A/«A/(0) ^—3 - e - b + 4e-^ -e" 26 + 6e-^ + ••• . (53) 

2y/nN2 1 J 

This result shows that, relative to the zero temperature value whose leading order is 

y/N (see (14T|) ). the temperature dependence of the net force Af starts at the order of 

1/N 2 . As a result, a noticeable deviation from A/(0) can be expected only for t ^> 1. 
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A/ fermions A / bosons 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 5. Low temperature approximations for A/ for fermions (left) and bosons 
(right), TV = 100. Thick solid line: the numerical result. Dotted line: the fermionic 
expansion (|53|) up to the term e~?, respectively the bosonic expansion (|58|) up to the 
term e _b . Dashed line: terms up to e~~, resp. e~ 2b . Dash-dotted line: terms up to 
e 2 ; resp. e . 



4-2. Bosons at low temperature 

For the bosons rj = 1, we know that at t — 0, the distribution is given by No = N and 
iVfc = for k > 0, and hence N = l/e a — 1 = N. This implies 

a = ln(l + i)«i (54) 

The net force at t = is thus 

A/(0) = A/ iVo = iV. (55) 

When the temperature grows from zero, a is expected to change continuously with 
t, and we may define the low temperature regime for bosons by a < 1 so that e~ a ~ 1 
along with t < 1. Then we can write for k > 0, 

iV fc = — i = 6 - w e~ bk + e~ 2bk + e~ 3bk + ■ • • , (56) 

where one can expect that, when n is high enough so that e~ na ceases to be near 1, the 
factor e~ nhk appearing in the nth term becomes very small. Consequently, in this low 
temperature regime, we have 

oo 

N = N - N k w N - e~ b - 2e~ 2b - 2e~ 3b . (57) 

k=l 

Plugging (|56j) and (!57|) into (13T|) and arranging terms according to the powers of e~ b , 
we find 

Af*N- 1 -e~ b - 9 -e- 2b - 1 ^e-* b + .... (58) 

We can see that, in the bosonic case, the temperature dependence of Af begins at the 
order of 1/N with respect to the zero temperature value. 
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5. The medium-temperature regime 

We have so far gained reasonably good approximations of the net force for high and low- 
temperature regimes. In this section we wish to study the temperature regime between 
the two, hoping to find some approximation to interpolate the previous two. Again, we 
consider the cases of fermions and bosons, separately. 



5.1. Fermions at medium temperature 

We have seen for fermions that, at t ~ 1, the net force A/(i) is still very close to A/(0). 
For a noticeable departure from A/(0), one needs to go above t ~ 1, presumably to 
t ^> 1 (or b <C 1), in which case the sum ([9]) can be approximated by the integral, 

e a+bk _ 1 Ay e a+ y* + 1 Ay Jo e a +y + 1 b ^ > v ; 

where we have used := bk, Ay := y k+1 — y k = b. This implies 

e~ a wet-1, (60) 

which is just ( I38jl . showing that the formula remains valid for medium temperatures as 
well. 

Precise evaluation of e~ a is particularly important in this temperature regime, 
because the final outcome of the net force Af is extremely sensitive to the variation of 
e~ a . For this reason, it is worthwhile to consider an improved approximation obtained 
from the fact that, for 'well-behaved' functions g(y) in [yo, oo) with ]hn y - i . 0O g(y) = 0, 
the trapezoid approximation of integrals yields 

Efld/^^ + T- /"%(y)dy, (6i) 
k=o 1 J y° 

which is better than that acting in ( 1591) by one order. As a result, one finds that the 
first correction term on the rhs of flBTj) improves (1601) to 

N-Nn/2 

e~ a w - 1. (62) 

This formula (|62|) reduces to (|6"0~1) for high temperatures where iVo <C 1 < iV, and 
to (1491) for low temperatures where N ^ 1 and the second term in (1621) is negligible 
compared to the first term. Note that the presence of N Q < 1 in the improvement is not 

N-l/2 JV 

insignificant, since at low temperatures e * can be significantly different from e~ on 
account of the ratio being 96 1. It also suggests that the deviation of iVo from its 
zero temperature value 1 may become important for medium temperatures, too. 

To acquire a meaningful formula between the low and high-temperature regimes, 
we can use fl62l) where the presence of N Q is expected to provide sensitivity to the low 
temperature regime, with the approximation that A^ is given by its high temperature 
value from ([60]), N ^ 1 - e~ N/t , that is, 

e (N-i^)/t _ L (63) 
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As can be seen in Figure EK this approximation formula for e~ a holds actually very well 
on the whole temperature regime. 

Similar improvement for the net force is a bit harder to achieve. Here, we will 
restrict ourselves only to leading-order approximations, and proceed as 

1 1 1 [b™ Ay 



fc =o ^ (k + 1/4) e &+bk - 1 Ay V 7T Jy k + 6/4 (e a +^ + 1) 
1 f 00 dy 



nb Jo J y _|_ 5/4 ( e <5+?, + i) 



(64) 



using fl30|) again, and considering 6 small enough. Introducing z := yy + b/4, we can 
rewrite it as 

A/-^=f - ** % 2 ^ d, 



For negative a, the resulting integral can be approximated by an asymptotic series 
[171 US] to obtain, 



Af « —=y/—a 



7TC 



[ 7T 2 1 77T 4 1 Z2. 
24 a 2 384 a 4 U e 



(66) 



The divergence of the formula for a — > may be dealt with by adopting the Pade form, 

A ^7^ 1 + ^ i 7 5? i > ( 67 ) 

which improves its validity towards a — > 0, 

In the medium-temperature regime, one may expect that the temperature is low 
enough for the second term of (1621) to be omitted, while it is high enough for N /2 to 
be omitted. If this is the case, one can simplify fl62|) to e~ a ps e~ or —a ~ y, which 
can be inserted into ( 1671) to obtain 

2 

1 + 2d f -t) + 23211 (_L) 
x ~ 24 \N J ~ 1152 V^/ 

Note that this formula takes care of the zero temperature value A/(0) correctly, and 
that, similarly to (1431) . the ratio Af/^/N is again a function of t and iV only through 
the combination t/N. 

Plotting the low+medium-temperature curve (|68|) . and the high temperature one 
( )43l) . we find that the two together actually cover the whole temperature range (see 
Figure [6b). From a technical viewpoint, we may find it convenient to introduce an 
interpolating function for the two curves connecting smoothly. This can be accomplished 
by using, for example, the scheme, 

/ \ 1- gi(x) + (x/x*) p ■ g 2 (x) 
flbtp(x) := ! + {xM p > ( 69 ) 

where g- m tp(x) gives the interpolating function of the two curves described by g%{x) and 
g2(x). The interpolating point may be chosen as the value where gi(a;*) = ^(^*) 
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Figure 6. Fermions, N = 100. a) The ratio of the approximate e~ a provided by 
(I6U)) to the numerically determined e~ a . The largest relative error is 10 -5 along the 
whole temperature region < t < oo. b) The net force Af. Solid line: the numerical 
result. Dashed line: the low+medium-temperature approximation (|68p . Dotted line: 
the high-temperature approximation Dash-dotted line: the interpolating curve 

with the choices mentioned in the text. 



holds, which ensures that at the weights 1+ ^ x p and j^r^y are equal irrespective 
of the value of p. 

In our case, we have x = t/N and choose the function Af in (1551) for gi and the 
function Af in (H3|) for g 2 . Then the numerically determined value of the interpolating 
point for N = 100 is found to be (t/N), = 0.237845. As for p, we know that flSHD is 
precise up to the order of (t/N) 5 , and if we maintain this then we find that the smallest 
choice is p = 5. Since this choice does not disturb the high-temperature expansion either, 
we may propose the resultant interpolating formula g intp in ( 1691) as the net force covering 
the whole temperature region. Figure [6b shows that this single formula reproduces the 
precise curve (obtained numerically) very well. Note that this interpolating procedure 
is actually independent of N, since, for all Ns, it involves the same functions (low and 
high temperature approximants of Af/y/N) of the identical variable t/N. 

5.2. Bosons at medium temperature 

This time we are again allowed to replace the sums by integrals with the approximation 
(loTj) . Since for low temperatures N is much larger than higher iV^, we preserve its 
discrete value and introduce the continuous variable only above k — 1. This gives 

iVi 1 r°° dy 



N = N + J2N k ^N + ^- + — f 
,-, 2 Ay Jv 



k=l 

1 h I 



AyJ yi e & +y-l 



e --i + ^i + 6 ln r^^- (70) 

For the medium temperature regime, we find it reasonable to assume that it is 
characterized by a C 1 and t ^> 1. Under this, one may simplify the outcome by 
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putting a = | in all terms except the first in the expansion (1701 (the advantage of the 
specific value \ will turn out soon), since by far the first term Nq = -^—r is the most 



4 „^ '"""yi ""—v- ^ Uii v> eQI _ 1 

sensitive term for the change of a on the scale of b or below. In this approximation, 
(T70|) leads to 



iVn 



AT 



i.n- 

6 1 



1 



6b 

4 



, 4 
tln-t. 

5 



e Q -l e x-l (v i-e 

The net force can also be evaluated by treating the sum analogously as 



(71) 



A/ = iV 



£A/ fc JV* 

k=l 



N + 



+ 



N + 
2 



A/iiVi 



dz 



- 1 



dz 



(72) 



2 e 4 — 1 V^b v , 

where we have gone through steps similar to ( 1651 and introduced an asymptotic 
expansion for the integrand to improve it in the dominant region. Note that the choice 
a — | has made this integral (^-independent. 

Actually, the approximation of the integrand with the asymptotic expansion 
becomes better if we keep only the part of the integral where the approximated integrand 
is positive. Because of this, we use 



rV2 


■ i r 


1 /~3> 

V 4 


.z* ~ 2. 



dz 



and (I7T|) to evaluate (1721 and find 

2 



In - + - 

5 5 5 



t - tint + 




(73) 



(74) 



We may further consider some corrections for the expression to render it a little nicer 
without sacrificing its precision. This is done by omitting the last term a/2 which is 
negligible compared to the first term N, which makes the t — limiting value the exact 
value A/(0) = N. We also replace the coefficient of t which is numerically 1.032 by 1 
for simplicity. We then end up with 



Af ^ N + t 




(75) 



As for lower temperatures the a-dependence of the terms k > 1 is suppressed more, 
and since we could reach a formula that is precise even at t = 0, we may hope that, 
in spite of the assumption a ~ | made above, this formula (1751 can be used even for 
the beginning part of the net force A/(t) including t = 0. Figure [7J shows that this is 
indeed the case. 

It is again possible to provide an interpolating formula between the low-temperature 
approximation and the high-temperature one based on the formula (1691) by choosing the 
functions (1751 and (|4"0~]) for gi and g 2 , respectively. Note, however, that this time the 
two curves do not cross each other, and the interpolating point £* may be chosen, e.g., 
as the location where the two functions differ the least. For example, for A^ = 100, 
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A/ 




t 



20 40 60 



Figure 7. The net force A/ for bosons, N — 100. Solid line: the numerical result. 
Dashed line: the low+medium-temperature approximation (|75[) . Dotted line: the high- 
temperature approximation (|40[) . Dash-dotted line: the interpolating curve (|69[) . for 
t* = 22.925 and p = 9. 

this yields i* := 22.925. In addition, the value of p can be chosen, for instance, by the 
requirement that the derivative of the interpolating function at t — be the same as 
that of the straight line that is the common tangent of the two curves to be connected. 
At N = 100, that common tangent straight line touches the low temperature curve at 
the point [7.338, 88.46] and the high temperature curve at [26.12, 45.66] in the coordinate 
plane of t and Af. The slope of that straight line is —2.278, and this is to be put equal 
to the derivative of flSS) at x m , i.e., % E (x*) = + J ^(*»)- 3 i(x*) ? where we have 

used = TpO^*) which comes from the definition of being the location where 

the difference is minimal. From this we find p = 8.641. One can invent some other 
criteria as well, but our aim here is only to present a formula in which the transition 
between the two curves appears as smoothly as possible, yet sharply enough to keep 
both curves practically intact on the regions where they are supposed to be reliable. 

Since the accuracy of the interpolation employed above is not very sensitive to the 
actual value of p, we may round it off to the nearest integer p = 9, say, for brevity. 
Figure [7] shows the resultant interpolating formula (1691) with our choice of functions, 
which is almost indistinguishable from the numerical result. We should, however, keep 
in mind that for bosons the interpolation must be done iV-dependently, because the 
scaling behaviour = 4^ ( j^j seen at the high-temperature regime does not arise 
at the low temperature regime. It also follows that a successful interpolation at some 
iV does not necessarily ensure a success for an analogously carried out interpolation at 
another N. 
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6. Conclusion and discussions 

In this paper, we studied the system of N particles confined in each of the two half 
lines separated by a partition wall at the centre of the harmonic oscillator potential. 
The partition is assumed to impose a set of distinct - the Dirichlet and the Neumann - 
boundary conditions on the left and on the right, respectively. Due to the discordance 
in the energy levels in two sides of the partition, and also to the different distributions 
of particles on the energy levels at finite temperatures, an induced force emerges on the 
partition, as we have seen earlier on the partition in the infinite potential well. 

We have evaluated the (dimensionless) net force A/(£) that arises on the wall 
both analytically and numerically as a function of (dimensionless) temperature t, and 
found that it exhibits a number of interesting behaviours characteristic to the harmonic 
potential. For instance, it has a non-vanishing limit A/(0) at the zero temperature 
limit, which is just N for bosons, while it is of the order of y/N for fermions. Note that 
in the case of the infinite potential well [9J, the value A/(0) for fermions is proportional 
to N 2 . 

On the other hand, in the high-temperature regime, the force A/(£) scales linearly 
in N for both fermions and bosons. As temperature grows, it decreases according to 
\j\ft and eventually vanishes in the limit t — ► oo. This is in sharp contrast to the 
infinite potential well case where the force A/(i) diverges according to \ft. 

The medium-temperature regime is somewhat difficult to deal with, but we have 
succeeded to obtain, after a rather technical argument, an analytic approximation 
that accounts for the numerical results reasonably well for both fermions and bosons. 
Interpolation to the low and high-temperature regimes can also be possible, and we 
presented a possible formula of the force covering the entire regime of temperature. 
Unlike the potential well case, the net force admits no minimum in the medium- 
temperature regime. 

The characteristic scaling behaviours in the zero temperature limit t — > can be 
understood heuristically. Namely, in the bosonic case, the force A/(0) comes entirely 
from the contribution of the ground level where all particles reside, and hence it is given 
by A/o multiplied by the number N of the particles. In the fermionic case, on the other 
hand, the force A/(0) consists of the contributions of Afk ~ (nk)"^ up to the Fermi 
level, yielding A/(0) oc J2k=i (x N^. This is to be compared to the infinite potential 
well case where we have A/* ~ k and hence A/(0) oc J2k=i k ^ N 2 . 

In the high-temperature regime, it can also be argued that the steady decrease in 
the present harmonic potential case, rather than the increase to infinity observed in the 
potential well case, derives basically from the spectral structure of the harmonic system. 
That is, for higher n the energy level difference between the two half harmonic systems 
remains constant and does not give larger contributions, in contrast to the potential 
well case where the energy level difference becomes larger and eventually diverges for 
n — > oo. In this respect, we recall the upper bound (|34l) for the force A/(£) whose 
decrease for higher t is already expected there as well. In physical terms, this is also 
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understood from the infinite stretch of the harmonic potential where the higher energy 
states can spread more in space. As a result, the energy of the system becomes less 
sensitive to a shift of the partition at the centre, resulting in the decrease in the force 
for high temperatures. 

What can we learn from the results obtained here for the harmonic system when 
combined with those obtained previously for the infinite potential well system? Suppose 
that a partition wall can actually be manufactured with the distinct set of boundary 
conditions assumed in this paper, and further that the induced force on the partition 
wall can be measured with sufficient accuracy. Then, one can estimate the profile of the 
potential in the neighbourhood of the partition wall by looking at the low temperature 
behaviour of A/, since the net force is sensitive to the spectral structure up to the Fermi 
level for fermions when the temperature is low. From the high temperature behaviour of 
Af, one can also obtain a crude picture of how the potential stretches in space further 
away from the wall. In addition, the characteristic scaling in the number N of particles 
in the zero temperature limit will reveal the statistics of the particles contained around 
the wall. In short, such a nontrivial partition wall may be quite useful in probing the 
profile of the potential as well as the statistics of the particles involved. The present 
probe is still primitive to use for generic potentials, but it can be improved if we learn 
further the behaviours of the net force for other types of potentials in the high and the 
low temperature regimes along with the scaling property in the particle number. 
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